arXiv: 1506.0838Ivl [quant-ph] 28 Jun 2015 


Reliability of optical gate array for C-Sign operator 

M.Lemus* Y.Ozhigovj N.Skovoroda^ 

June 30, 2015 


03.65, 87.10 


Abstract 

We study optical gate array of KLM model of quantum computer for C- 
Sign operator, which contains linear elements and two nonlinear phase shifts. 
Linear elements and photon counting are taken to be ideal, whereas nonlinear 
phase shifts are subject of realistic and unavoidable factors of decoherence: 
the immediate photon leak from cavities and mediated leakage with the atom, 
which flies out of the cavity and plays the role of ancilla. We use JC model 
with Lindblad master equation to estimate the dependence of the validity of 
C-Sign operator from the rate of photon leak and estimate the limitation, 
which this validity imposes to quantum computations. 


1 Introduction and background 

Quantum mechanics gives us the new paradigm of computations: the peculiar type 
of information processing realized by quantum computers, called quantum compu¬ 
tations. It has been shown by [T] and [2] that for some specihc tasks, quantum 
computation allows for faster performance than it is theoretically possible in a clas¬ 
sical computer. Quantum information and encryption provides a secure way of 
transmitting information based on the wavefunction collapse property of quantum 
systems and the use of quantum encryption algorithms. One of the main challenges 
of the held is the physical realization of a quantum computer. One of the proposals, 
the KLM scheme 0010. based on linear optical components, has paved the way 
for different implementations based on quantum optics. 

The photon, being a chargeless and massless particle, is a very adequate candi¬ 
date as the carrier of quantum information. The advantages of using photons include 
an easy and low-loss way to transport them using optical hber, their relatively well 
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understood quantum properties, and the existence of well-developed techniques for 
generating and measuring photons. 

Linear quantum optics has proven to be very good playground for experimental 
and theoretical developments in the realization of quantum computers. The pho¬ 
tons allow for at least two viable representations of quantum information in form 
of the dual rail and the polarization representation. The parametric down conver¬ 
sion [6] and the single atom microlaser [7j allow us to produce single photons on 
demand, ready for computation. In addition to that, the KLM scheme provides a 
mechanism from which all the single qubit gates can be constructed using only mir¬ 
rors, phase shifters and beamsplitters, which are all common and well understood. 
Even the nonlinear effects needed to build the CNOT gate can be obtained by the 
non-deterministic nonlinear shift (NS) gate. The associated decoherence of those 
components is low, meaning that most of the information loss comes from the pho¬ 
ton storage between operations and from the photon measurement. Both of these 
can be reduced using optical resonators to store light in the infrared and microwave 
domain, and using solid state photoelectric detectors based in semiconductors to 
measure the qubits. 

The main problem with the KLM scheme is scalability; this means that the re¬ 
sources needed for computation should increase linearly with the number of qubits 
used. The fact that the NS gate has a low success probability makes successive appli¬ 
cations of it unreliable, as the probability of success quickly decreases exponentially. 
This effect can be countered using quantum teleportation gates as described in [8], 
but the amount of entangled pairs of qubits needed to make the CNOT gate work 
with 95% reliability is in the order of 100 pairs per gate |9] which limits the scheme 
scalability. 

Another solution to the scalability problem is to extend the scheme to include 
nonlinear components. One way of realizing the NS gate is using nonlinear Kerr 
media to produce cross phase modulation, theoretically solving the main problem of 
the optical quantum computer model. Nevertheless, there is always some absorption 
associated with this nonlinearity, and it has been estimated HD] that in the best 
known case, approximately 50 photons must be absorbed to get one to experience 
a TT cross phase modulation. This means that currently, the use of Kerr media isn’t 
a solution for the scalability problem. 

Cavity QED poses itself as a solution to the challenge of scalable qubit entan¬ 
glement, both for qubits represented as atom states with interaction mediated by 
photons as shown by HU, as well as for qubits represented by photons interacting 
through atoms inside a cavity. A single atom in a cavity is enough to produce an 
entangled state between photons, as shown experimentally by Recently, a setup 
by [l3] realizes the NS gate with an upper bound error probability of approximately 
2.47%. This work continues by generalizing it to explore how the inclusion of decou¬ 
pling and dissipative effects change the potential reliability of the NS and C-Sign 
gates. 
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Figure 1: NFS: optical cavity, which the two-level atom flies through. The energy 
of the held in the cavity does not exceed 2uc- We choose the appropriate time 
tq = As/v for hnding the atom in the cavity for the realization of NFS. 

2 Description of the physical system 

Consider an optical Fabry-Perot cavity (as shown in Figure 1) with a two-level 
Rydberg atom hying through it with speed v. The geometrical parameter w denotes 
the half-width of the beam waist at the center of the cavity and L denotes the 
distance between the mirrors. Let the low and high energy states of the atom be 
denoted by \g) and |e) respectively an the n photon state inside the cavity be \n). 
Suppose that the atoms inside the cavity are prepared initially in the \g) state, 
and the system is built so no more than two photons can enter the cavity. The 
Hilbert space associated with the system is S'panjl^f, 0), l^f, 1), l^f, 2), |e, 0), |e, 1)}. 
The evolution of such system is described by the Jaynes Cummings Hamiltonian 


(H): 


Hjc = oJcO^a + uJaOzI^. 'y{a~^a + cxa^) 


( 1 ) 


where af,ai are operators of creation and annihilation of photons; a^,a are 
operators of excitation and relaxation of the atom in the cavity, and 7 is the coupling 
constant for interaction between the atom and the held {h is taken to be 1). This 
Hamiltonian leaves invariant subspaces of the form A„ = Span{\g,n + 1), |e, n)} and 
the space = Span{\g, 0)}. By taking the detuning A = Uc — tOa, Hjc, limited on 
A„ can be written as a 2 x 2 matrix 



1 0\ / A /2 

0 1 )^ -^/2 


Hjc = {n + l/ 2 )a;c 
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The evolution of the system can be fully obtained by solving the eigenvalue 
problem for Hjc, which gives the eigenvalues for the energy: 

e± = {n + l/2)u, ± nj2 (3) 

where + A'y'^{n + 1) is the generalized Rabi frequency of the system. 

The associated eigenstates are the dressed states of the system 

|+,n) = cos{9n)\g,n + l) + sin{9n)\e,n) 

|-,n) =-sm( 6 '„)| 5 (,n + 1 ) + cos( 6 '„)|e,n) 

where tan{9n) = {2\/n + l 7 )/(f 2 „ — A). For the system ground state 1 ( 7 , 0 ), the 
energy eigenvalue is given by 


eg = (-l/2)n;e-A/2 (5) 

Considering the system as a part of an array of components, in which the number 
of photons is hxed, we turn our attention to only the interaction and free atom part 
of the Hamiltonian. It is useful then, to consider the interaction picture associated 
with Hq = {a'^a + l/2)a;c. In such picture, the evolution of the state according to 

0 is 


i){t) = aoe + ^(a+je j) + a.je j)) (6) 

j 

Let the system be initially prepared in a state of the form -0(0) = ao|( 7 , 0 ) + 
ai\g, 1) + a 2 \g, 2). Then, equation 0 shows that after a time t the system evolves 
to the state 


V’(t) =aoe**^/2|^,0) + 

+ ai[{cos{yLQt/2) — icos{29Q)sin{yLQt/2))\g, 1 ) — 

— ism(26*o)sm(f2oV2)]|e, 0)+ (7) 

+ a2[{cos{Vtit/2) — icos{29i)sin{VLit/2))\g, 2) — 

— isin{29i)sin{VLit/2y\\e, 1) 

The dissipation inside the cavity can be handled using a master equation in the 
Lindblad form on the system’s density matrix p 

ihp =[H,p]+iY^ i^j{LjpL+ - ^{L+Lj, p}) (8) 

j 

The main source of decoherence inside an optical cavity is the photon leakage 
through the mirrors, the effect of which can be described by the photon annihilation 
operator acting on the electric held state Li = a and the photon escape rate constant 
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Figure 2: C-Sign gate array 




Kphi which depends on the Q factor of the cavity and the wavelength of the photons. 
The resulting Lindblad equation is 

ihp = [Hjc, p] + iKphiapa'^ - ^{a+a, p}). (9) 

3 Mathematical model of C-Sign gate array 

The C-Sign gate array, proposed in [13] is shown in Figure 2. The wires are denoted 
by Xi,X 2 ,yi, 2/2 and each wire can contain 0,1 or 2 photons of the fixed frequency ojc 
and circular polarization, with a single limitation: the number of photons that are 
present in the gate array simultaneously does not exceed 2. The number of photons 
in any wire we denote by the same letter. 

Decoherence comes from the non ideal work of photo detectors, for example, 
dark countings, photon leakage from the walls of the cavities and from the harmful 
entanglement between photons and atom, which can remain when the atom flies 
out of the cavity. In this work we ignore the errors of photo detectors and analyze 
only the photon leakage and harmful entanglement. We also assume that the linear 
optical elements work ideally, thus the single source of errors is connected with the 
work of optical cavities. 

The logical states of our computations are the spatial photon states in the dual 
rail representation. The qubits are denoted by x and y in Figure 2 and the compu¬ 
tational basis is defined as 


|0) = |0>,|1)2 |l> = |l)i|0)2 (10) 

The atomic states play the role of ancilla, which is prepared in the state \g) and 
in the ideal case must leave the gate array in the same state. Unfortunately, the 
ideal situation is unattainable in the given model even theoretically. We can only get 
an increasingly accurate approximation to it by choosing larger times of interaction 
between the field and the atom in an optical cavity. However, this increase of time 
will inevitably lead also to the increase in the probability of photon leakage out of 
the cavity. Our task is to find the optimal time for finding the atom in an optical 
cavity, depending on the intensity of photons leak. 

The array consists in four wires, two idle ones X 2 ,y 2 and two wires xi,yi which 
are acted upon by the optical components. The idle wires cannot be disregarded, 
as the photons inside of them contribute to the overall energy of the system. 
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There are two linear gates: beam splitters BSl and BS2, whose action in Fock 
space looks as 


h)ai|m)ai = 


\/n\m\ 


{atn4r\^)aM 


Ci2 


+ 4)I”[4w - aj)l”|0>..|0>„, 

yn\m\ yl yl 


( 11 ) 


and two optical cavities, which perform the nonlinear shifts: NSl, NS2. The 
nonlinear shift in the ideal case acts on Fock states of photons as 


| 0 )^| 0 ), | 1 >^| 1 ), | 2 )^-| 2 > ( 12 ) 

This is realized by Jaynes-Cummins Hamiltonian Q. The gates work sequen¬ 
tially, from the left to the right, as shown in the Figure 2. We assume that the beam 
splitters work instantly, whereas the NS gates work simultaneously, and their work 
lasts for the time At. Nonlinear shifts are realized by optical cavities with resonant 
frequency cUc, which is close to the frequency of atom transition from the excited to 
the ground state uja with a small detuning A. The atom, initially prepared in the 
ground state, flies through the cavity along the path, and during the flight interacts 
with held inside the cavity via ([^. We interrupt this interaction by the dehation of 
photons from the cavity in the moment t = t that is chosen in advance; the choice 
of r determines the accuracy of the gate array. 

At this stage we meet the unavoidable decoherence factor, which comes from the 
residual entanglement between atom and held, mentioned above. 

To choose r we consider three arithmetic progressions 


A . fO fO 

JS. . 6 q , 6]^, . . . , 

R • 

JD . 6 q , 6^, . . . , 

O . Lq, 6;^, ..., 


ti- 

'-''m 1 • 




(13) 


where A, B and C consist of the sequential time moments, when the unitary evo¬ 
lution of JC transforms photonic states as |0) to |0), |1) to |1) and |2) to —12) 
respectively. If we choose the time r, when these progressions approximately coin¬ 
cide, we will obtain the required duration. 

The evolution of Fock states is shown in ([^. Eigenvalues of the Hamiltonian 
blocks (|^ for n = 0,1, 2 have the form 

e, = (-l/2)a;e - A/2 

e±,o = (l/2)a;, ± V'A^T^/2 (14) 

e±,i = (3/2)a;, ± + 87 V 2 

Nonlinear dependence of amplitude from the photon number, called optical non¬ 
linearity is a rare and exceptional phenomenon. A manifestation of this exclusivity 
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is the fact that the match in progressions (13) is only approximate. We could con¬ 
sider the generalization of operator ( 12 ): | 0 ) —)■ e“| 0 ), | 1 ) —)■ e*^|l), | 2 ) —)■ e*‘^| 2 ) 

, where a -|- c 7 ^ 26, and realize its particular case with only one separated optical 
cavity with flying atom by the choice of appropriate parameters A and the time r. 

The nonlinearity means that for some f > 0, integers kg^k±fl,k±^i and a -|- c 7 ^ 
2b (modulo 27r) (nonlinearity of the phase on photon number) the following four 
equations are true: egt = 2TTkg + a, e±fit = 27rfc±^o + b, e±^it = 2TTk±^i + c. 

By inspection of equation ([^, it is clear that and must be real numbers, 
which means that 6 , c G {0, tt} (modulo 27r). We note that for the zero detuning (A = 
0) the impossibility of the exact coincidence between elements of the sequence B and 
C is evident, because \/2 is not rational and thus the differences of these progressions 
are incommensurable. We can only choose the value of r that is maximally close to 
the both progressions. The more accuracy we wish to obtain the larger r will be. 
We are only interested in the accuracy of the scheme, not the time of its execution. 
The value r can thus only affect the quality of the circuit via the second factor of 
decoherence - the leak of photons through the walls of cavities. 

For the detuned case, the nonlinearity condition can be satished if we take such 
d = A/7 so that \/4 -|- d‘^/\/8 + is rational; that can be guaranteed if we hnd 
the rational solutions of the equation = 1 , for which there is dense set of 

appropriate values for the parameter d. To achieve nonlinearity it is then sufficient 
that a ^ c (modulo 27r). That we can achieved by almost any arbitrary choice of 
parameters among these solutions. For the C-Sign case we need more, (12) is realized 
only when a = c + it (modulo 27r). The numerical computations with Mathematica 
shows that there is no solution of such a system. The best we can do is to hnd the 
approximate solution. 

The evolution with the leak of photons gives us the mixed state even if the 
initial state was pure; we thus must use the density matrix p for the description of 
the evolution. We use master equation 

The evolution of the state in the work of our gate array thus can be described 
by the equation ([^, where the Hamiltonian H and Kph-, Hat are piecewise constants; 
namely, H coincide with the Hamiltonian, corresponding to the action of linear 
elements (for the simplicity, we can apply the corresponding transformation 0 
directly in the corresponding time) without dissipation, and in the time frame when 
the held interacts with the atoms we set H = Hjc and take the actual value of the 
K coefficients. 

We use the distance S(po,p) = ||po ~p|| between the ideal density matrix po and 
the real density matrix p at the end of computation as the criterion of the validity 
of the gate array. This value d depends on the intensity k of photon leak, and we 
can hnd this dependence for the diherent choice of tq. 
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4 Numerical simulation of the density matrix dy¬ 
namics 

All gates, except for the main part of the NS gates, are considered being ideal (both 
beamsplitter gates and the phase shift inside the NS gates, if enabled). 

The NS gates simulation was performed in the assumption that both of them 
have identical time. The simulation operates with 6 particles (4 rails and 2 atoms): 
Xii X 2 , Hi, 1 / 2 , of, a2. There are 19 total possible states in the system, because all 
the other states are unreachable when using a valid input. The simulation deals 
only with those 19 reachable states. 

4.1 The Hamiltonian for the NS gate 

The Hamiltonian has to take into an account all the particles in the system. x2 and 
y2 receive only the phase shift over time with an amplitude which is assumed to be 
equal to Uc- Detuning A is equal to Ua — ooc- 

The total Hamiltonian for the main part of the NS gate (excluding the phase 
shift) has the following form; 

Hns= Ucixtx];+xtx2 +ytyi+yty2)+^aiaia];+ata2)+ 

g{afx]^+xfai+a2X2+X2a2) ^ 

Variable g could be excluded from the parameters list, if all the variables are 
taken relative to g. When doing so, the main parameters for the calculation are the 
following: Wa/g,A/g (or A/wa), t = T\/2g/7r (where T is the absolute NS gate 
duration), h is taken equal to 1. 

Also one extra parameter is phs, that denotes enabling [phs = 1) or disabling 
{phs = 0) a phase shift gate at the end of both NS gates. 

From here on, it is taken g = 0.1 (does not affect anything as all parameters are 
relative to g), wc = g * * 10®. Values for detuning A that are used below are 

from 0 to 5g, which is hve orders of magnitude lower than the cavity frequency. 

4.2 Lindblad operators for photon leakage 

We simulate the photon leakage from the cavities using the Lindblad master equa¬ 
tion. With the inclusion of the Lindblad master equation in the diagonal form, the 
simulation equation takes the following form; 

= UstPtUst + — \{L*LiP + pL*Li)), 

i ■ 6t 




( 16 ) 



Lindblad operators for the photon leakage inside the NS gates have the following 
form: 

L\ ^ ly ^ 37 ^ 

^2 ~ ^ ^2 

In the following equations ly is the photon leakage coefficient. For the simulation, 
we take it relative to the photon-atom interaction strength g, so this adds one more 
parameter to our model: ly/g. 



4.3 Calculating the error rate 

Error rate error takes values from 0 to 1 and is equal to 1 minus the gate validity. 
The gate validity from here on is calculated as max{\eigenvalues{Pexpected — Presuit))) 
where Ptest = Kl^i) + l^/i) + + \y 2 )){{xi\ + {yi\ + (xsl ( 2 / 2 I) density matrix was 

used as the input. Also other random valid Prand density matrices with the total 
number of photons equal to 2 (gate restriction) were tested as input, and the result 
corresponded with the result for Ptest- 


4.4 The case of zero detuning 

In the case of zero detuning, the only viable parameter sets are those with whole 
values of t. That is explained in section and could be also seen from the hgure 



Figure 3: Error over the selected NS gate time, with and without the phase shift, 
without detuning. Optimal values are 3, 7, 17, 41, 99. Photon leakage is disabled. 
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The higher the NS gate duration is, the more precise match we can hnd, and 
the lower is the error rate, but as could be seen on hgure the faster does the error 
rate grow over increasing the photon leakage and/or the inaccurate detuning setting 
(which is the same as the detuning itself for the case when the desired detuning is 
equal to 0). There are no reasons for choosing higher durations over lower ones, 
if the gate validity does not improve, so we treat as “optimal” only those sets of 
parameters that lower the error rate compared to the best found parameters set 
with lower time values. 



Error over detuning Error over photon leakage 

(logarithmic scale on leakage) 

Figure 4: Error rates over detuning and photon leakage for the optimal cases. 

4.5 The case of non-zero desired detuning 

It was also shown that non-zero detuning levels could greatly improve the validity 
of the gate. On hgure it could be seen how does the error rate depend on the 
absolute detuning value. The graph is done only for low times for clarity, the picture 
is the same for higher times, but the range of detuning values that interests us grows 
with time and also becomes more hne-grained. 
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Figure 5: Error rate over detuning and time. 

It is visible that for values of t around 4, some non-zero detuning range actually 
lowers the overall error level. In the case of non-zero detuning, t is not required to 
be a whole number, and it is not whole for most of the optimal parameters set. 

If parameter sets with non-zero detuning value is taken into consideration, it 
greatly widens our possible selection of the “optimal “ parameter sets, and also im¬ 
proves the lowest reachable error rate for a given time range. This could be seen on 
figure Also the same figure shows the overall dependency of the reachable error 
rate over our duration range. 
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Figure 6: Optimal time values and their error rates for the cases with and without 
the desired detuning. 


Figure shows how does the error rate in the cases of zero detuning depend on 
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the inaccurate detuning setting A — Aopt and from the photon emission coefficient. 
It could be seen that the dependency is not stronger than the corresponding depen¬ 
dency on hgure so non-zero detuning does not introduce larger error values in 
cases of inaccurate detuning setting and/or non-zero photon leakage from the NS 
gates cavities. 

With t = 99, the lowest reachable error rate is ~ 0.008 without detuning and 
0.002 with A 4.80 * g. 



error 



Error over detuning offset 


Error over photon leakage 
(logarithmic scale on leakage) 


Figure 7: Error rates over detuning and photon leakage for the optimal cases with 
desired detuning. 


5 Conclusions 

We have conducted an estimation of the error rate of modihed KLM model. Optimal 
time values were found for the zero detuning case, and it was shown how non-zero 
detuning could lower the error rate and broaden our choice for the time values. 

Error introduced by the photon leakage from the cavity at the NS gates was 
simulated, and it was shown how the photon emission rate affects the optimal time 
value. 

Error introduced by inaccurate detuning setting was simulated, it was shown 
that it depends mostly on the chosen time. 

Optimal {time, detuning) pairs were found for the cases of non-zero detuning, 
and it was shown that those are more efficient (introduce a lower error for the same 
or lower duration) compared to the cases when the detuning is equal to zero. 

It was also shown that the error introduced by photon leakage from the cavities 
and by the inaccurate detuning setting imprecision does not increase in the non¬ 
zero optimal detuning case when compared to zero detuning case. Thus, to get 
the optimal gate validity for the known photon escape rate and detuning setting 
imprecision, one could hrst calculate the maximum NS gate duration after which 
the error introduced by them raises to an unacceptable level, and then hud the 
optimal values for t and ^ where t is slightly less than that critical value. 
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The choice of the duration tq can be made by exact hxation of the moment to 
when atom entries the cavity; for the speed of atoms about 1cm. per sec. taken from 
a reservoir of cold gas such hxation can be performed with high precision. 

Graphs at the hgures iHI shows that the lowest relative error level accessible 
by the considered model of C-Sign is about 0.01; for the increasing noise it rapidly 
increases. If we roughly estimate for what number of qubits Grover search algorithm 
([2]) can be realized on this model, taking into account the linear growth of error 
in course of quantum computation and take the error on one Grover operator equal 
to G-Sign error we found, we obtain the number n ^ 13. This estimation, coming 
for the account of different sources of decoherence is more pessimistic than the 
more abstract estimations (see, for example, [H]) that gives the values n ^ 20. To 
reach such values in the realistic models, like KLM, we need the following additional 
modernization of this scheme. 
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